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We analyze the entanglement between two modes of a free Dirac field as seen by two relatively 
accelerated parties. The entanglement is degraded by the Unruh effect and asymptotically reaches 
a non- vanishing minimum value in the infinite acceleration limit. This means that the state always 
remains entangled to a degree and can be used in quantum information tasks, such as teleportation, 
between parties in relative uniform acceleration. We analyze our results from the point of view 
afforded by the phenomenon of entanglement sharing and in terms of recent results in the area of 
multi-qubit complementarity. 
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I. INTRODUCTION 

Entanglement plays a central role in quantum informa- 
tion theory. It is considered a resource for quantum com- 
munication and teleportation, as well as for various com- 
putational tasks 0. The importance of understanding 
entanglement in a relativistic setting has received con- 
siderable attention recently 0, 0, S 13 ■ Such an under- 
standing is certainly relevant from a fundamental point 
of view, since relativity is an indispensable component 
of any complete theoretical model. However, it is also 
important in a number of practical situations, for exam- 
ple, when considering the implementation of quantum 
information processing tasks performed by observers in 
arbitrary relative motion. 

Entanglement was shown to be an invariant quantity 
for observers in uniform relative motion in the sense that, 
although different inertial observers may see these corre- 
lations distributed among several degrees of freedom in 
different ways, the total amount of entanglement is the 
same in all inertial frames In non-inertial frames, 
entanglement was first studied indirectly by investigat- 
ing the fidelity of teleportation between two parties in 
relative uniform acceleration More recently, the 

observer-dependent character of entanglement was ex- 
plicitly demonstrated by studying the entanglement be- 
tween two modes of a free scalar field as viewed by two 
relatively accelerated observers 

A uniformly accelerated observer is unable to access 
information about the whole of spacetime since, from 
his perspective, a communication horizon appears. This 
can result in a loss of information and a correspond- 
ing degradation of entanglement. In essence, the accel- 
eration of the observer effects a kind of "enviromental 
decoherence" , limiting the fidelity of certain quantum 
information-theoretic processes. A quantitative under- 
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standing of such degradation in non-inertial frames is 
therefore required if one wants to discuss the implemen- 
tation of certain quantum information processing tasks 
between accelerated partners. 

In curved spacetime two nearby inertial observers are 
relatively accelerated due to the geodesic deviation equa- 
tion. Accordingly, the results of indicate that in 
curved spacetime even two inertial observers will disagree 
on the degree of entanglement in a given bipartite quan- 
tum state of some quantum field. Indeed, a thorough 
investigation into entanglement in an expanding curved 
spacetime shows that entanglement can encode informa- 
tion concerning the underlying spacetime structure |^. 

In this paper we analyze the entanglement between 
two modes of a Dirac field described by relatively ac- 
celerated parties in a flat spacetime. We are interested 
in understanding how both the crucial sign change in 
Fermi-Dirac versus Bose-Einstein distributions, and the 
finite number of allowed states in fermionic systems due 
to the Pauli exclusion principle (in contrast to the un- 
bounded excitations that can occur in bosonic systems), 
affect the degradation of entanglement produced by the 
Unruh effect. We find that unlike the bosonic case, where 
the entanglement degrades completely in the infinite ac- 
celeration limit, in the fermionic case the entanglement 
is never completely destroyed. We analyze the degra- 
dation of entanglement in the system by applying the 
constraints of entanglement sharing |^ and track the in- 
formation originally encoded in these quantum correla- 
tions by employing a set of multi-qubit complementarity 
relations 0. As in Q, our results can be apphed to the 
case that Alice falls into a black hole whilst Rob barely 
escapes through eternal uniform acceleration. 

The remainder of this paper is organized as follows. 
In Sec. ^ we consider two modes of a free Dirac field 
that are maximally entangled from an inertial perspec- 
tive. Two parties, an inertial observer named Alice and 
a uniformly accelerating observer named Rob, are each 
assumed to possess a detector sensitive only to one of the 
two modes. Each measures the field with his/her detec- 
tor and the results are compared in order to estimate the 
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entanglement between the modes. 

Section Hill discusses the Unruh effect for Dirac parti- 
cles as experienced by Rob. If a given Dirac mode is in 
the vacuum state from an inertial perspective, then Rob's 
detector perceives a Fermi-Dirac distribution of particles. 
This has a strong effect on the entanglement that exists 
between Alice and Rob, and therefore plays an impor- 
tant role in any quantum information task they might 
perform that uses this entanglement as a resource. 

In Sec. II VI we calculate the entanglement between the 
modes from the perspectives of both Alice and Rob. Due 
to the presence of a Rindler horizon, Rob is forced to 
trace over a causally disconnected region of spacetime 
that he cannot access. Accordingly, his description of 
the system takes the form of a two-qubit mixed state. 
We calculate the entanglement using mixed state entan- 
glement measures such as the entanglement of formation 
[3 and the logarithmic negativity ,9,]. We also estimate 
the total correlatio ns ( classical plus quantum) via the 
mutual information Our results show that the en- 

tanglement of formation does not vanish as it does in the 
bosonic case, but rather reaches a minimum of I/V2 m 
the limit that Rob moves with infinite acceleration. 

Since the fermionic system we are considering is accu- 
rately described by a pure state of three qubits, we study 
the constraints placed on the system by the phenomenon 
of entanglement sharing in Sec. Our analysis shows 
that no inherently three-body correlations are generated 
in the quantum state. That is, all of the entanglement 
in the system is in the form of bipartite correlations, re- 
gardless of Rob's rate of acceleration. 

Using complementarity relations applicable to an over- 
all pure state of three qubits, as well as to the various 
two-qubit marginals, we identify the different types of 
information encoded in the quantum state of our system 
in Sec. IVII This enables us to study how specific subsys- 
tem properties depend on Rob's rate of acceleration and 
to explain how some of the entanglement from the iner- 
tial frame is able to survive in the fermionic system, even 
at infinite acceleration. Finally, we summarize our re- 
sults and suggest possible directions for further research 
in Sec. EUl 



II. THE SETTING 

Consider a free Minkowski Dirac field in 3 + 1 dimen- 
sions 

17'" 9^1/; — mip ~ 0, 

where m is the particle mass, the Dirac gamma ma- 
trices, and ^ is a spinor wavefunction. Minkowski coor- 
dinates = (ct,x) with fi — {0,1,2,3} are the most 
suitable to describe the field from an inertial perspec- 
tive. The field can be expanded in terms of the positive 
(fermions) and negative (antifermions) energy solutions 
of the Dirac equation ip'^ and tp'^ respectively, since they 



form a complete orthonormal set of modes, 

j dk{au^t + bli,^). (1) 

In the above. A: is a notational shorthand for the wavevec- 
tor k which labels the modes. The positive and negative 
energy Minkowski modes have the form 
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where ujk — {m? + k^)^/^, and tps = 0s (k) is a constant 
spinor with s — {tii} indicating spin-up or spin-down 
along the quantization axis, satisfying the normalization 
relations [TJ ±0^ = {ujk/m) 5ss' , 4>t' 4>ti = 0, with 
the adjoint spinor given by = (jif^ 7°. The above posi- 
tive and negative energy solutions satisfy the orthogonor- 
mality relations 

{^t ,^t') = ~ i^k ) - '5(fc - fc') , ii^k M') = 0, 

where the Dirac inner product for two mode functions is 
given by 

{(f){x,t),Lp{x,t)) = I dx(t>\yL,t)ip{y.,t). 



The modes i/'^ are classified as positive and negative fre- 
quency with respect to (the future-directed Minkowski 
KiUing vector) dt for ujk > 0, i.e. 



k ' 



UJk > 0. 



The operators aj, , 6^ and , are the creation and an- 
nihilation operators for the positive and negative energy 
solutions of momentum k which satisfy the anticonmuta- 
tion relations 

{ai,a]} = {bi,b]} = dij, 

with all other anticommutators vanishing. The 
Minkowski vacuum state is defined by the absence of any 
mode excitations in an inertial frame. 



io)=niO'' 
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where the — } superscript on the kets is used to indi- 
cate the particle and anti-particle vacua, respectively so 
that afelOfe)"*" = bk\Ok)^ = 0. We will use the notation 
here, and throughout the rest of the work, that the mode 
index (k,k' . . .) will be a subscript affixed to the occu- 
pation number inside the ket, and that the absence of a 
subscript on the outside of the ket indicates a Minkowski 
Fock state. Since (a^)^ = (&jj)^ ~ 0, there are only two 
allowed states for each mode, |0a:)+ and |1a:)+ = aj,|0fc) + 
for particles, and similarly for anti-particles. 

Consider two maximally entangled fermionic modes in 
an inertial frame, 

\h.,k,) = ^(|OfcJ + |O.J+ + |lfcJ + |Uj + ), (2) 
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where the subscripts A and R indicate the modes as- 
sociated with the observers Ahce and Rob, respec- 
tively. AU other modes of the field are in the vac- 
uum state, and therefore the state can be written as 
I*) = \'l^k.,kn)mk^k,,kjOk)+Uk'\Ok')-}- Now as- 
sume that Alice is stationary and has a detector sensitive 
only to mode fc^. Rob moves with uniform acceleration 
and takes with him a detector that only detects particles 
corresponding to mode kn. We ask the question of what 
is the entanglement between modes kj^ and kji observed 
by Alice and Rob, given that Rob undergoes uniform ac- 
celeration. Note that in order to determine the amount 
of entanglement, Alice and Rob perform measurements 
which are then compared by either party in order to es- 
timate the correlations in the results. Due to Rob's ac- 
celeration, at some point Alice's signals will no longer 
reach Rob, but Rob's signals will always be available to 
Alice, (see Fig.QJ. At this point only Alice can compare 
the measurement results and estimate the entanglement 
of the state. Let us now consider the state observed by 
Rob. 



III. UNRUH EFFECT FOR DIRAC PARTICLES 

Consider Rob to be uniformly accelerated in the {t, z) 
plane (c = 1). Rindler coordinates (r, C) are appropriate 
for describing the viewpoint of an observer moving with 
uniform acceleration. Two different sets of Rindler coor- 
dinates, which differ from each other by an overall change 
in sign, are necessary for covering Minkowski space 
These sets of coordinates define two Rindler regions that 
are causally disconnected from each other. 

at = e"^sinh(ar), az = e"^ cosh(ar), in region I 

(3) 

at = — sinh(aT), = — e"^ cosh(ar), in region II 

where a denotes Rob's proper acceleration. The above set 
of coordinates both give rise to the same Rindler metric 

ds^ = dt^ - dz^ - d^Xi = e^"? (dV - d'^C) - d^^^, 

where xj^ = {x,y) are the same in both Minkowski and 
Rindler spacetimes. A particle undergoing eternal uni- 
form acceleration remains constrained to either Rindler 
region I or II and has no access to the opposite region, 
since these two regions are causally disconnected. Fig- 
ure serves to illustrate these ideas, as well as to intro- 
duce our labeling scheme, where we refer to the acceler- 
ated observer in region I as Rob (R) and to the corre- 
sponding fictitious observer confined to region II (whose 
coordinates are the negative of Rob's) as anti-Rob (R) . 

The coordinates (t, () have the ranges — oo < t, C < oo 
separately in region I and in region II. This implies that 
region I and II each admit a separate quantization proce- 
dure with corresponding positive and negative energy so- 
lutions {V'fc'*'; V'fc"} and {'>Pk^'^ Since the Rindler 
metric is static (independent of r) it will admit solutions 
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FIG. 1; Rindler spacetime diagram: Lines of constant po- 
sition are liyperbolas and lines of constant proper time r 
for the accelerated observer run through the origin. Note 
that while r flows in the direction of t in region I, it flows in 
the direction of ~t in region II, (i.e the dashed line rotates 
counter-clockwise for increasing values of r). A uniformly ac- 
celerated observer Rob (R) with acceleration a travels on a 
hyperbola constrained to region I, while a fictitious observer 
anti-Rob (R) travels on a corresponding hyperbola in region 
II given by the negative of Rob's coordinates. The horizons 
H± are lines of r = ±cxd which Alice (A) will cross at finite 
Minkowski times. 

of the form ej*"'^ (/)q(C,xx), with ipa a spatially depen- 
dent spinor . Particles and anti-particles will be clas- 
sified with respect to the future-directed timelike Killing 
vector in each region. In region I this is given by dr where 

dt dz 
dr = T-dt + ^ a {zdt + tdz), 

OT OT 

which is a boost into the instantaneous comoving frame 
of Rob. Thus, mode solutions in region I having time 
dependence ipl'^ ~ e"*"^"^ with w > represent positive 
frequency solutions since dripi~^ = However, 
in region II dr points in the opposite direction of dt (in- 
creasing r flows in the direction of —t, see Fig.nj. Hence 
in region II the future-directed timelike Killing vector is 
given by d-r — —dr 0,0. Thus, a solution in region II 
with time dependence e""^''" with w > is actually a neg- 
ative frequency mode since d-r e~"^'^ = iliJ e""^''". Hence, 
the positive frequency mode in region II is given by 
^]^+ ^ e^^^ with > satisfying d-ripk^ = -iw<pl^+ . 
Due to the causally disjoint nature of region I and II, 
the modes V'^^ have support only in region I and van- 
ish in region II, while the opposite is true for the modes 
-0^^^ in region II. The Rindler modes satisfy orthonor- 
mality relations analogous to the Minkowski modes [l^ 

= 0- and (Vfe*,Vfe'^) = ±5,,,,5{k-k') 
where a S {/, //}. 

In region I, let us denote (c^,c^^) as the annihila- 
tion and creation operators for fermions (particles) and 
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as the annihilation and creation operators for 
anti-fermions (anti-particles). The corresponding parti- 
cle and anti-particle operators in region II are denoted 



as {cl' , I. J. 

anti-commutation relations {c^,c^/^} = {d%,d'J^i^} — 



and (dj. 



These obey the usual Dirac 

Skk' , with all other anti-commutators, including 
those between operators in region I and the causally dis- 
connected region II, equaling zero. Taking into account 
the two sets of modes in each Rindler region, the Dirac 
field can be expanded, in analogy to Eq. (^), as 



dk 



(4) 

Eq. 1^ and Eq. Q) represent the decomposition of the 
Dirac field in Minkowski and Rindler modes, respectively. 
We can therefore relate the Minkowski and Rindler cre- 
ation and annihilation operators by taking appropriate 
inner products. Using the Rindler orthogonality rela- 
tions and Eq. I^J we have — (V'fe , 4') ■ Substituting 
Eq. for ijj in this last expression yields 



= / dk 



' ere {/,//} (5) 



where the Bogoliubov coefficients are given by the in- 
ner product of the Rindler modes wavefunctions with the 
Minkowski positive and negative frequency modes 



(6) 



A similar calculation for yields the corresponding ex- 
pression 



dk 



'kk' 



kk' 4' 



a e {1,11} (7) 



with the same Bogoliubov coefficient as in Eq. In 
deriving Eq. |(7J) use has been made of the following prop- 
erties of the Dirac inner product: {(j)i,(j)2)* = (0*, 02 ) = 
{4>2, 4>i)- The calculation of the Bogoliubov coefficients is 
straightforward, though lengthy and an exercise in spe- 
cial functions. Details can be found elsewhere 23. EM 13 ■ 
For our purposes, the end result of such calculations 
yields a relationship between the Minkowski and Rindler 
creation and annihilation operators given by the Bogoli- 
ubov transformation 



ak 




. blk . 





-e sinr 
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(8) 



where tanr = exp(— 7r57) with Q = uj/{a/c), the ratio of 
the frequency lu to the only naturally occurring frequency 
in the problem a/c, and is an unimportant phase that 
can always be absorbed into the definition of the opera- 
tors. It is easy to see from Eq. (0 and its adjoint that 
given the anti-commutation relations of the Rindler op- 
erators, the Minkowski anti-commutation relations are 
preserved. In Eq. ||SJ we have made the single mode ap- 
proximation 0, which is valid if we consider Rob's de- 
tector as sensitive to a single particle mode in region I 



such that we can approximate the frequency uja observed 
by Alice to be the same as the frequency as observed 
by Rob, i.e lja ^ ujr = to Note that this Bogoli- 

ubov transformation mixes a particle in region I and an 
anti-particle in region II. Correspondingly, the Bogoli- 
ubov transformation that mixes an anti-particle mode in 
region I and a particle in region II is given by 



e ^ smr 
-e smr cosr 



cos r 





\ < 1 







(9) 



Since the anti-commutators between particle and anti- 
particle operators and between region I and region II 
Rindler operators are zero, it is easy to see that the 
Minkowski operators in Eq. (jSjl anti-commute with the 
Minkowski operators in Eq. @, as they should (since k 
and —k represent two separate modes). 

As stated above, Eq. ((Sj) reveals that the Minkowski 
particle annihilation operator is a Bogoliubov transfor- 
mation between a particle in region I of momentum k 
and an anti-particle in region II of momentum — fc, i.e. a 
transformation that mixes creation and annihilation op- 
erators. We can understand this in terms of our previous 
discussion of the time dependence of positive frequency 
Rindler modes in regions I and II. For a massless Dirac 
field a positive frequency Rindler mode has the form 191 
~ exp(iA:(^ — iwr) (for the scalar case see 0, Ir^ l 
and a positive frequency Rindler mode in region II has the 
form i/'fc^^ ~ exp(ifcC -|- iwr). Thus, in order to construct 
a positive frequency Minkowski mode ?/'^ that extends 
analytically from region I to region II, we need a lin- 
ear combination of t/;^^ and (V'i^fc^)* so that both have 
the space and time dependence exp(jfc^ — iioT). This 
means that the Minkowski operator must mix and 

4]_l. as in Eq. JHJl. Similarly, to construct a Minkowski 
mode that analytically extends '^]}^ from region II to 
region I we must form a linear combination of "^^^^ and 
k^^~k)* both have the space and time dependence 

exp(i/cC -|- iuir). This corresponds to 6^ mixing d\. and 
cl'l as in Eq. For a massive Dirac field the time 

dependence of the Rindler modes remains the same, but 
the spatial dependence is more complicated '^l- The 
above argument for the mixing of Rindler particles (anti- 
particles) in region I (II) and anti-particles (particles) in 
region II (I) goes through unchanged. 

Having related Minkowski and Rindler creation and 
annihilation operators, we now wish to relate the 
Minkowski vacuum to the corresponding Rindler vacuum. 
It is useful to note that Eq. © can be written as a two- 
mode squeezing transformation |2l| 
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S 



4^1 
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for the single mode k with S given by 



exp 



d^Ile 
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Using the relation B e"-^ = B+[B,A] + [B, [B,A]]/2\+ 
[B,[B,[B,A]]]/3l + ■■■ and the identity [AB,C] ^ 
A {B, C} - {A, C} B we find 



acting on the lower expression for \0k)'^ in Eq. H15() . yields 
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2 /t 

cos r Cj, ' 



and 



SciS^ 



cos r Ci. 



cl 



re 
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k)ll 



= [cos^ r 4^+ Sin' rciULid'Jl) |0,); |0 



e-''^sinrd"^ 



(10) 



= cj\0kn\0-k)n, 



(17) 



where in the second equality we have used the anti- 
commutation relations between cj} and d^Jf, to obtain 



cos r d^Jj^ 



e '"^ sin rc{. 



(11) f 



Now afc and 6_fe, respectively, annihilate the sin- 
gle mode particle and anti-particle Minkowski vacua 
OfclOfc)^ = and 6_fc|0_fc)~ = 0. Since mixes particles 
in region I and anti-particles in region II, we postulate 
that the Minkowski particle vacuum for mode k in terms 
of Rindler Fock states is given by 



|Ofe)+ = ^A„|nfe)+|n_fe),-, 



where '22 



ci^\Okn = \ik)r 



-k 



|0- 



(12) 



(13) 
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= li- 



fe///- 



As a comment on notation, the Rindler region I and II 
Fock states carry a subscript I or II, respectively, on the 
kets while the Minkowski Fock states are indicated by 
the absence of an outside subscript on the kets. Momen- 
tum mode labels are attached to the Fock occupation 
number, and the — } ket superscript indicates a par- 
ticle or anti-particle state, respectively. Applying ak from 
Eq. Iini) to Eq. lO, we have 



= flfelOfe)" 



cos r Cl. 



e sin r d^^l 



1 

E 

n=0 



An \nk)] \n-k) 



II 



{Ai cosr-^oe"'* sinr) |Ofe)+ |l_fc)7^ 
Ai = Aoe""* tanr. 



(14) 



Normalization +(0fc|0fe)+ [Aop- 
cos r so that we finally arrive at 



I All 2 = 1 yields Aq 



|Ofe)+ = cosr|0fe)+|0_fe)7/ 
COST + e""'' 



-^^ d"^ 



sinr|U)+|l_fe),-, 



where the second equality in Eq. (|15() is very useful for 
keeping track of any possible transposition signs aris- 
ing from the anti-commuting operators when applying 
Minkowski operators to |Ofe)'^. As an important exam- 
ple that we will use subsequently, a short calculation of 
al |Ofe)+ using the adjoint of Eq. ((117)) 



COST c 



s+^'^sinrd- 



// 
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a minus sign upon transposition, and in the third equal- 
ity we have used the anti-commutation relations to write 



1 — d^J^ d[^j, , noting that the latter term an- 
nihilates the anti-particle vacuum. Lastly, one can easily 
verify that a\. acting on |1a;)^ yields zero. 



.|lfe>+ = (c 



H<ti 



sin r d^Ji 



ensuring that (a^)^ 



:) \lk)i\Q-k)u=0, 

(18) 



0. 



We note that the form of the Dirac particle vacuum 
for mode k in Eq. ((T5|l . which can be written as |Ofc)+ = 

cos r J2n=o t^-ii" ^ Wk)i \n-k)Jj, is complementary to the 
form of the Minkowski charged scalar vacuum for mode 
k, which is given by 0,^3 



= (coshr)"^ ^ tanh" 



r \nk) 



II 



where in the latter case r is defined by tanhr = 
exp(— Trfi). Qualitatively, in going from the scalar field 
to the Dirac field, scalar mode functions are replaced by 
spinors, the infinite number of equally spaced bosonic 
levels are replaced by two fermionic levels, and the hy- 
perbolic functions in the Bogoluibov transformation be- 
tween Minkowski and Rindler modes/operators are re- 
placed by the corresponding trigonometric functions (in 
essence r — > i r in going from the scalar to the Dirac 
field). 

The two Minkowski states, lOfc)"*" and |lfc)^, corre- 
spond to the particle field of mode k observed by Alice. 
On the other hand, an observer moving with uniform ac- 
celeration a in one of the regions has no access to field 
modes in the causally disconnected region. Therefore, the 
observer must trace over the inaccessible region, consti- 
tuting an unavoidable loss of information about the state, 
which essentially results in the detection of a mixed state. 
Thus, when a Minkowski observer detects a vacuum state 
|0fc)"'"(0fe| for mode fc, an accelerated observer in region I 
sees a distribution of particles according to the marginal 
state describing region I, 



pi 



= Trii [|0fc) + (0fc|] 
= cos2r|0fc)++(0fc 



• 2 I 

sm r 



Ik) I I 



(19) 



As the region I observer accelerates through the 
Minkowski particle vacuum |Ofe)"*" of mode k his detec- 
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tor registers a number of particles given by 



(Ofe|c^^c^|0,)+ = rr,, 



= Tri 



Pk 



i^ciiOknOk 



sin^r J{lkKci\lk)r 



1 



sm r = 



1 



(20) 



where use has been made of tanr = exp(— Trfi) with Q = 
ujc/a, and we have defined the Unruh temperature (where 
ks is Boltzmann's constant) as 



T 



ha 
ks^Trc 



(21) 



Equation H20(l is known as the Unruh effect which 
shows that the uniformly accelerated observer in region 1 
detects a thermal Fermi-Dirac (FD) distribution of par- 
ticles as he traverses the Minkowski vacuum. Qualita- 
tively, we can understand the Unruh effect as follows. 
The constant proper force F that acts on a mass m 
(a detector) in Rob's instantaneous co-moving frame, 
to keep it under uniform acceleration, can be written 
as = ma (which can be integrated as dp/dr — ma 
where p = jmv, 7 = (1 — f^/c^)~^/^, to yield the hy- 
perbolic orbits in Eq. ©). The work 5W that is per- 
formed on the particle to keep it in uniform accelera- 
tion is the product of F times a characteristic distance 
through which the force acts, which we take to be the 
Compton wavelength of the particle 6z — h/{mc). This 
yields SW — FSz = (ma) h/{mc) = ha/c, which is pro- 
portional to ksT in Eq. (|21|1 . Thus, the energy that 
is supplied to keep m under constant acceleration goes 
into exciting Rob's detector, and curiously has a ther- 
mal spectrum (see 2'||). In the general relativistic case 
this energy is supplied by the gravitational field acting 
on a stationary observer (Rob) situated outside a black 
hole, who experiences constant acceleration by nature of 
his stationarity. A freely falling observer (Alice) who 
eventually crosses the event horizon would experience the 
Kruskal vacuum (the Minkowski-like vacuum of Fig. ^ 
and thus detect no Unruh radiation with a co-moving 
detector. 

Note that in the case of a scalar field, Eq. (|20() has a 
crucial minus sign in the denominator leading to a Bose- 
Einstcin (BE) distribution of particles appropriate for 
bosons. In both the scalar and Dirac case, the state of the 
Minkowski vacuum for mode fc is a two-mode squeezed 
state (bosonic and fermionic, respectively), which to a 
uniformly accelerated observer confined entirely to re- 
gion I, is detected as a thermal state (BE and FD, re- 
spectively). 

In the following we investigate how the Unruh effect for 
Dirac particles affects the entanglement between various 
Dirac modes. Specifically, we study the entangled Bell 
state given by Eq. ^ in the case that Rob is uniformly 
accelerated in region I. Assuming a detector for Rob that 
is sensitive to a single Rindler particle mode (fc ku), 
we decompose Rob's single particle Minkowski states in 



Eq. (0) into the appropriate Rindler particle and anti- 
particle states utiHzing Eqs. l(T^ and ((T7|l . We then pro- 
ceed to evaluate various measures of entanglement. As 
stated in the introduction, the advantage of utilizing a 
Dirac field over a bosonic scalar field is that due to the 
finite occupation of the fermionic states, we obtain finite 
dimensional density matrices that lead to closed form 
expressions for the entanglement measures that are more 
easily interpreted than their infinite dimensional bosonic 
counterparts. 



IV. FERMIONIC ENTANGLEMENT FROM A 
NON-INERTIAL PERSPECTIVE 



In this and subsequent sections we will use the follow- 
ing notation: A will indicate the inertial observer Alice, 
I will indicate the uniformly accelerated observer Rob 
(R) confined to region I, and II will indicate the ficti- 
tious complementary observer anti-Rob (R) in region II 
(see Fig. Pi, that arises from the second (negative) set of 
Rindler coordinates in Eq. ||3J). Furthermore, since we are 
in the single mode approximation, we will drop all labels 
(fc, — fc, fcyi, . . .) on states and density matrices indicating 
the specific mode. Thus, the Minkowski particle mode 
|nfc)+ for Alice will be written \n)A, the Rindler region 

I particle mode \nk)^ |n)/, and the Rindler region 

II anti-particle mode \n-k)n !"•)// (for n G {0,1}). 
Likewise, we will refer to the "Minkowski mode for Al- 
ice" simply as "mode ^" , the Rindler particle mode in 
region I as "mode I," and the Rindler anti-particle mode 
in region II as "mode II." 

The density matrix for the Minkowski entangled state 
in Eq. (01 is, from an inertial perspective 



inertial 

Pa,r 



/ 1 1 





V 1 1 



in the basis |00), |01), |10), |11) with \ab) = |a)A |6)_r. 
Here, and only here, we have used the subscript R to 
indicate Rob's Minkowski Fock states in Eq. 

To describe the entanglement of the state as seen by 
an inertial Alice and a uniformly accelerated Rob, we 
expand the Minkowski particle states \0)r and |l)i^ into 
Rindler region I and II particle and anti-particle states 
using Eqs. ifT^ and ((T7|l to obtain 



1 



PAJJI 



'r|000)(000| + sin2r|011)(011| + |110)(110|) 



= - (COS 

+ i (cosrsinr|000)(011| -|-cosr|000)(110| 
+ sinr|011)(110| +h.c.). 
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which we write in matrix form as 







cos^ r 








cos r sin r 








cosr 


0\ 
























































1 




cos r sin r 








• 2 

sm r 








sinr 





2 




























































cos r 








sinr 








1 







I 























0/ 



PA,I = 77 










cosr 



• 2 

sm r 



cos r \ 





1 / 



(22) 



in the basis |00), |01), |10), |11) where \ab) = \a)A \b)i. 

To determine whether or not this state is entangled 
we use the partial transpose criterion . This criterion 
specifies a necessary and sufficient condition for the ex- 
istence of entanglement in a mixed state of two qubits. 
If at least one eigenvalue of the partial transpose of the 
density matrix is negative, then the density matrix is 
entangled. The partial transpose is obtained by inter- 
changing Alice's qubits {\aAbi){cAdi\ — > 6/)(aA 



T 

Paj 



( 



V 





sm r cos r (J 

cos r 

1 y 



J has eigenvalues (1, 1, cos^ r, — cos^ r)/2, the last of 
which A_ = — i cos^ r, is always negative for < a < cx) 
i.e. for < r < 7r/4. This means that the state is always 
entangled. 

Quantifying entanglement for mixed states is fairly in- 
volved (2^]. A pure state of a bipartite system a and 
h can always be written in the Schmidt basis, |f/')ab = 
^„ a„|n)o|n)f,, where the quantum correlations between 
the states are evident. Entanglement between the sys- 
tems is given by the von-Neumann entropy of the reduced 
density matrix pa defined as S{pa) — —tr{palog2{pa)), 
which is a function of the Schmidt coefficients a„. Un- 
fortunately, there is no analog of the Schmidt decompo- 
sition for mixed states, and the von-Neumann entropy is 
no longer a good measure of mixed state entanglement. 

A set of conditions that mixed state entanglement mea- 
sures should satisfy is well known [2^ . There is no unique 
measure, and several different mixed state entanglement 
measures have been proposed. Among the most popu- 
lar are those related to the formation and distillation of 



entangled states. Consider the number m of maximally 
entangled pairs needed to create n arbitrarily good copies 
of an arbitrary pure state using only local operations and 
classical communication. The entanglement of formation 
is defined as the asymptotic conversion ratio, m/n in the 
limit of infinitely many copies @, 



EpiPab) = min^KS'(p^), 



in the basis |000), |001), |010), |011), |100), |101), |110), |111), 
where for notational convenience we have defined 
\abc) = \a)A\b)i\c)ii. Note that the unimportant phase 
factor (j> discussed in Section Hill has been absorbed into 
the definition of the creation and annihilation operators. 

Since Rob is causally disconnected from region II, we 
take the trace over the mode in this region, which results 
in a mixed density matrix between Alice and Rob 



where the minimum is taken over all the possible realiza- 
tions of the state pab = EiKl*ab)(*a&l- 

The opposite process gives rise to the definition of the 
entanglement of distillation. This is the asymptotic rate 
of converting non-maximally entangled states into max- 
imally entangled states by means of a purification pro- 
cedure. The entanglement of distillation is in general 
smaller than that of formation. This shows that the en- 
tanglement conversion is irreversible and is due to a loss 
of classical information about the decomposition of the 
state. Bound entangled states are a consequence of this: 
no entanglement can be distilled from them even though 
they are inseparable. 

The entanglement of formation can be explicitly cal- 
culated for two qubits. It is given by 



Ef= ~ log2 



1 - VI - C2 , M - VI - C2 
^ log2 n 



where C is the concurrence. The concurrence of a pure 
state lijj) of two qubits is given by 

where ay (Jy IV'*) represents the 'spin-flip' of \ip) and 
"*" denotes complex conjugation in the standard basis. 

The generalization of the concurrence to a mixed state 
p of two qubits follows by minimizing the average con- 
currence over all possible pure state ensemble decomposi- 
tions of p, defined by a convex combination of pure states 
Si = {pi,tpi}, such that p = J2iPi\'^i){'^i\- ^^^^ ^^Y^ 



i 



Wootters |3] succeeded in deriving an analytic solution 
to this difficult minimization procedure in terms of the 
quantities defined as the square roots of the eigen- 
values (which are all positive) of the non-Hermitian op- 
erator pp, ordered in decreasing order. The matrix 
p — (ay (8) cTy) p* {ay ® a y) Is the spin- flip of the quan- 
tum state p which not only exchanges the states |0) and 
|1), but, in general, also introduces a relative phase. The 
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closed form solution for the concurrence of a mixed state 
of two qubits is given by 

C (p) = max {0, Ai - A2 - A3 - A4} , Xi > Xi+i > 0. 

Unfortunately, the entanglement of distillation cannot be 
explicitly calculated in this case. Therefore we will use 
the logarithmic negativity which serves as an upper 
bound on the entanglement of distillation. It is not an 
entanglement measure because it does not satisfy the re- 
quirement of being equal to the von-Neuman entropy for 
pure states. However, it is an entanglement monotone 
since it satisfies all other criteria to quantify entangle- 
ment. It is defined as N{p) = log2 ||p"^||i where ||p"^||i 
is the trace-norm of the partial transpos e density matrix 
i.e., the sum of the eigenvalues of yvP^V^- Since 
the matrices p^ in this work are symmetric, N{p) is sim- 
ply given by log2 of the sum of the absolute values of the 
eigenvalues of p^ . 

To quantify the entanglement of pAj in Eq. H22|l we 
compute the spin-flip matrix pAj 



PAJ = 7; 



( 



sin^ 



cosr \ 





y cos r cos^ r J 



and find that 



cos-^ 



PAJ PAJ = 2 






cos r 








cos^; 



has eigenvalues (cos^ r, 0, 0, 0). The concurrence is then 
given by C{pAj) = Ai = cosr, which is unity at zero 
acceleration, as expected, and approaches the value 1/V2 
for infinite acceleration r — > 7r/4. The entanglement of 
formation is 

„ 1 . N , . 1 + sin r , 
Ef= - -(l+smr)log2( ) 

. s, 1-sinr 
- -(1 -smr)log2( ), 

and the logarithmic negativity is 

N = log2(l + cos^ r) 

In an inertial frame, r = 0, and the state of the sys- 
tem defined by Eq. ||2Jl is maximally entangled. In the 
limit of infinite acceleration, the logarithmic negativity 
is log2 I = 0.585, implying that the entanglement in the 
infinite acceleration limit is finite. This means that the 
state is always entangled and can be used as a resource 
for performing certain quantum information processing 
tasks. This is in strong contrast to the bosonic case P|, 
where it is found that the entanglement goes to zero in 
this limit. As in |3| the infinite acceleration limit can 



be interpreted as Alice falling into a black hole while 
Rob hovers outside the black hole, barely escaping the 
fall. Since close to a black hole horizon spacetime is flat, 
Rob must be uniformly accelerated in order to escape 
the black hole. However, Alice falls in as she is an iner- 
tial observer. The connection to the present situation is 
made by recognizing that, according to Rob, a commu- 
nication horizon appears causing him to lose information 
about the state in the whole of spacetime. As a result, 
the inertial entanglement is degraded, and there will be 
a reduction in the fidelity of any information processing 
task performed by Alice and Rob using this state. The 
analogy with the black hole scenario is further strength- 
ened by observing that classical information can only fiow 
from Rob to Alice after Alice has crossed the horizon. 

We may also calculate the total correlations between 
any two subsytems of the overall system by using the 
mutual information 10|, 

I = S{pa) + S{pa) - S{pab)- 

This measure quantifies how much information two cor- 
related observers (one with access to subsystem a and the 
other with access to subsystem b) possess about one an- 
other's state. Equivalently, it represents the distance be- 
tween the actual joint distribution and the product state 
obtained when all correlations are neglected. Calculat- 
ing the relevant marginal density operators, the mutual 
information of the state H22|) is found to be 



1 „ 
1 = 1 — - cos r log2 



1 



(23) 



1 



(1 - - cos^ r) log2 2 

1 , , ^ , f 1 + cos^ r 
-(1+cos^r) log2 ( 



i(l -cosV) log2 



1 — cos r 



where we have used the fact that the density matrix in 
region I is p/ = i [cos^ r|0)/(0| + (sin^ r + so 
that 



0/ ^ 1 2 , cos^ r 

S{pi)= - 2 cos ?-log2(^— 



(24) 



(1 



logo 1 



2 ' °^ V 2 
Similarly, the density matrix of mode A is pA = 
i(|OA)(OA| + |1a)(1a|), yielding S{pa) = 1. Finally, us- 
ing Eq. |(22Jl we find that 

S{PAJ)= - i(l + cos^r)log2(^i±|^) (25) 

1 2 ^ 1 1 - COS^ ^ 

- ""^"^ ^°S2 ( 2 

The mutual information is equal to 2 at r = and goes 
to unity in the infinite acceleration limit. This behavior 
is reminiscent of that seen in the bosonic case j^]. 



9 




FIG. 2: (Color online) Bipartite pure state entanglement. 
Thick solid curve (Blue): Anti-Rob in region II with (Alice 
and Rob in region I). Thin dashed curve (Red): Alice with 
(regions I and II). Dotted curve (Green): Rob in region I with 
(Alice and anti-Rob region II.) 



B. Mixed state entanglement 

For mixed state entanglement we can consider, besides 
the entanglement between the Minkowski mode A and 
Rindler mode I calculated in Sec. IIVI two additional bi- 
partite divisions of the subsystems. 

• The entanglement between mode A and mode II. 
Tracing over the mode in region I, we obtain the 
density matrix 



PAji = 



( 





sin^ r 
sinr 






sinr 
1 









The partial transpose of pa,ii is given by 



V. ENTANGLEMENT IN OTHER PARTITIONS 
AND ENTANGLEMENT SHARING 

To explore entanglement in this system in more detail 
we consider the tripartite system consisting of the modes 
A, I, and II. In an inertial frame the system is bipartite, 
but from a non-inertial perspective an extra set of modes 
in region II becomes relevant. We therefore calculate 
the entanglement in all possible bipartite divisions of the 
system as well as any tripartite correlations that may 
exist. 



A. Pure state entanglement 



T _ ^ 

Pa,ii - 2 





sinr 



sinr ^ 

sin^r 

1 





which has eigenvalues (1,1, sin r, — sin r)/2, the 



last of which A_ = — i sm 



r, is less than or equal 
to zero. At r = the eigenvalue is zero, which 
means that there is no entanglement at this point. 
However, for any r > entanglement does exist 
between these two modes according to the partial 
transpose criterion. The logarithmic negativity in 
this case is given by = log2(l 



sm^ r ) 



Calculating the spin-flip of pA,ii 



We investigate pure state entanglement in three dif- 
ferent bipartite divisions of the system. Since the over- 
all state of the system is pure, these entanglements are 
uniquely quantified by the von-Neumann entropy. In 
what follows we make repeated use of the fact that, 
for a bipartite pure state pab with subsystems a and 6, 
S(pab) = and S{pa) = S{pb)- There are three separate 
cases to consider. 

• The entanglement between mode A and (modes I 
and II) is given by S{pa) = S{piji) = 1. Thus, 
these two subsystems are always maximally entan- 
gled, regardless of the value of the parameter r. 

• The entanglement between mode I and (modes A 
and II) is given by Eq. since S{pi) = S{pA,ii)- 

• Finally, the entanglement between mode II and 
(modes A and I) is given by Eq. H25|l since S{pii) = 

S{PAj)- 

Figure |(2J) shows the entanglement in each of these bi- 
partite partitions as functions of r. 



PAJI = 2 



we find that 



/O 




sin r sin 















cos^ r 



PAJI PAJI = 2 



/O 0^ 

sin^ r sin'^ r 

sin r sin^ r 

\ J 



has eigenvalues (sin^ r, 0, 0, 0). Thus, the concur- 
rence is given by C{pAji) — \i — sinr, which 
is zero at zero acceleration as expected, and ap- 
proaches the value 1/^/2 for infinite acceleration 
r — > 7r/4. The entanglement of formation is 



T-. 1/1 M /l + cosr 
Ef= - -(l + cosr)log2 ( 

- ^(l-cosr)log2 f — 1^ 
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and the mutual information is 

/ = S{pa) + S{pii) - S{pAji) 

1 / 9s, / 1 + cos^ r 
= 1- -(l + cos^r)log2 ( 

1 2 m 1 ~ cos^ r 

- -(1 - cos^r)log2 



icos2rlog2 



(1 - i cos^ r) log2 (l - ^ cos^ ^ 



At r — the mutual information is zero and ap- 
proaches unity as the rate of acceleration goes to 
infinity. 

The entanglement between mode I and mode II. 
Tracing over the modes in A, we obtain the density 
matrix 



which has eigenvalues (2 sin^ r, 2 cos^ r, 1 + (1 + 
sin^ 2r)i/^ 1 - (1 + sin^ 2r)i/2)/4^ the last of which 

A_ — — 1[1 — Vl + sii?2r] is less than or equal 
to zero. Again, at r = there is no entanglement 
between these two subsystems. Yet, similar to the 
last case, entanglement does exist between these 
two modes in non-inertial frames according to the 
partial transpose criterion. Further, the logarith- 
mic negativity N — log2(5(l + \/l + sin^ 2r)) is 
nonzero for all r > 0. The spin- flip of pi // is given 

by 



Piji = 2 



/ sin^ r cos r sin r \ 

10 



y cos r sin r cos^ r J 



and the matrix 



PIJI = 2 



^ cos^ r cos r sin r \ 








1 







'sinr sin r J 



PijiPiji = 2 



/ cos^ r sin^ r 



y cos r sin^ r 



cos r sin r 





cos^ r sin^ r J 



The partial transpose of piji (obtained by inter- 
changing /'s qubits) is given by 







^ cos^ r 








' ) 




T 


1 








cos r sin r 







PlJI = 


2 





cos r sin r 


1 











I 








sin^ r J 


1 



has eigenvalues (cos^ r sin^ r, 0, 0, 0). Thus, the 
concurrence is given by C{piji) = Ai = sinr cos r 
which is zero at zero acceleration, and approaches 
the value 1/2 for infinite acceleration r tt/A. The 
entanglement of formation in this case is 



Ef= - 



(1 + \/ 1 — sin^ r cos^ r 



) , / 1 + V 1 — sin'^ r cos^ r 
- log2 o 



(26) 



(1- VT 



• 2 9 

sm r cos^ r 



sm r cos^ r 



and the mutual information is 



S{pi) + S{pii) - S{piji) 

1 2 1 f cos^ r 
--cos r log2 ( — ^ 

(1 - i cos^ r) log2 (^1 - ^ cos^ r 

1 , 9s, / 1 + cos^ r 
-(1 + cos^r) log2 ( 

I . 2s, / 1 ~ cos^ r \ 
-(1-cos r) log2 ( 1-1. 



Again, we find that the mutual information is zero at 
r = 0, and increases to a finite value (in this case 
/ = 3/2 (2 — log2 3)), as the rate of acceleration goes to 
infinity. The results of the above calculations are shown 
in Figs. ® - ©. 



C. Tripartite entanglement 

Entanglement in triparite systems has been studied by 
Coffman, et. al., Q for the case of three qubits. They 
found that such quantum correlations cannot be arbitrar- 
ily distributed amongst the subsystems; the existence of 
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Negativity 




r 



8 



FIG. 3; (Color online) The logarithmic negativity as a func- 
tion of r. Thick solid curve (Blue): between Alice and Rob 
in region I. Thin dashed curve (Red): between the modes in 
regions I and II. Dotted curve (Green): between Alice and 
anti-Rob in region II. 



E_F 




"8" T 

FIG. 4: (Color online) Entanglement of formation as a func- 
tion of r. Thick solid curve (Blue): between Alice and Rob 
in region I. Thin dashed curve (Red): between the modes in 
regions I and II. Dotted curve (Green): between Alice and 
anti-Rob in region II. 



Mutual Information 
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FIG. 5: (Color online) Mutual information as a function of 
r. Thick solid curve (Blue): between Alice and Rob in region 
I. Thin dashed curve (Red): between the modes in regions I 
and II. Dotted curve (Green): between Alice and anti-Rob in 
region II. 



three-body correlations constrains the distribution of the 
bipartite entanglement which remains after tracing over 
any one of the qubits. For example, in a GHZ-state, 
I GHZ) = 1 000) -I- I 111), tracing over one qubit results in 
maximally mixed state containing no entanglement be- 
tween the remaining two qubits. In constrast, for a W- 
state, |W) = |001) + |010) + |100), the average remaining 
bipartite entanglement is maximal. Coffman, et. al., 
analyzed this phenomenon of entanglement sharing 
using an entanglement monotone known as the tangle, 
defined as the square of the concurrence t = . They 
also introduced a new quantity, known as the residual 
tangle, in order to quantify the irreducible tripartite cor- 
relations in a system of three qubits (a, &, and c) |^. The 
definition is motivated by the observation that the tangle 
of a with b plus the tangle of a with c cannot exceed the 
tangle of a with the joint subsystem 6c, i.e., 

Ta,b + 'ra.c<Ta(b^c)- (27) 

Subtracting the terms on the left hand side of Eq. H27(l 
from that on the right hand side yields a nonnegative 
quantity referred to as the residual tangle Ta^b,c, i-e., 

Ta,b,c = 'ra{b,c) — Ta,b — Ta.c- (28) 

The residual tangle or 'three tangle' is interpreted as 
quantifying the inherent tripartite entanglement present 
in a system of three qubits, i.e., the entanglement that 
cannot be accounted for in terms of the various bipar- 
tite tangles. This interpretation is given further support 
by the observation that the residual tangle is invariant 
under all possible permutations of the subsystem labels 

3- . . . 

To quantify tripartite entanglement in our system we 

use the residual tangle ta,i,ii = ta(^iji) — ta,i — ta,ii- 
This quantity is zero for the situation we are considering, 
since taj = cos^ r, ta,ii = sin^ r and ta(iji) = 2(1 — 
Tr[pj jj]) = 1. Thus, the state that we are studying has 
no tripartite correlations for any value of the acceleration 
rate. Instead, any entanglement existing in the system is 
necessarily bipartite in nature. 

The marginal bipartite tangles, plotted in Fig. as 
functions of r, are found to be strongly constrained. For 
low rates of acceleration, modes A and I remain almost 
maximally entangled while there is very little entangle- 
ment between modes I and II and between modes A and 
II. As the acceleration grows, the entanglement between 
modes I and II and between modes A and II increases, 
while the entanglement between modes A and I is de- 
graded. The main system of interest (mode A plus mode 
I) becomes increasingly entangled to mode II and there- 
fore, after tracing over mode II, we observe an effect anal- 
ogous to environmental decoherence in which the modes 
in region II play the role of the environment. The distri- 
bution of entanglement in the system is non-trivial since 
entanglement between modes A and I is not conserved as 
it is in the inertial case. 
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VI. COMPLEMENTARITY 

We next analyze our system in terms of several com- 
plementarity relations designed to identify the different 
types of information encoded in a quantum state in an at- 
tempt to better understand how various subsystem prop- 
erties depend on Rob's rate of acceleration. Specifically, 
we are interested in explaining (i) the nonconservative 
nature of the entanglement discussed in the last section 
and (ii) the fact that not all of the initial entanglement 
between Alice and Rob is destroyed, even at infinite ac- 
celeration. The latter of these two results corresponds to 
the most obvious difference between the fermionic case 
studied here and the bosonic case investigated in P|. 

We begin by making use of the relationship 



77(p) + r(p) + 52 {pa) + S^ {Pb) = l, 



(29) 



which shows that an arbitrary state of two qubits (a 
and b) exhibits a complementary tradeoff between the 
amounts of separable uncertainty 77, bipartite entangle- 
ment (as quantified by the tangle) r, and a unitarily in- 
variant measure of information about the single particle 
properties that it encodes. The separable uncertainty 



r; (p) = Tr (pp) + M {p) ^ t [p) , 



(30) 



5: vip) ^ 1, is a measure of the uncertainty or ig- 
norance encoded in the two-qubit mixed state p regard- 
ing individual subsystem properties that is unrelated to 
the presence of entanglement between the qubits. In 
Eq. (|30|l . p represents the spin- flip of p, and M {p) = 

1 - Tr is the marginal mixedness of p. 

Similarly, (pk) is a measure of the informa- 
tion pertaining to a single qubit encoded in the 
marginal density operator pk- Specifically, {pk) = 
1/2 \y'^ (pk) (pfe)] is the average of the squares of 
the single qubit properties associated with qubit k = a,b. 
The first of these properties, the coherence v of qubit k, 
quantifies, e.g., the fringe visibility in the context of a 
two-state system incident on an interferometer, and is 
given by 



where a 



(fc) 

+ 



1 




Tr 



(fc) 

Pk<y+ 



(31) 



is the raising operator acting on 



qubit k. Similarly, the predictability p which quantifies 
the a priori information regarding whether qubit k is in 
the state |0) or the state |1), e.g., whether it is more likely 
to take the upper or lower path in an interferometer, is 
given by 



p(Pfc) 



Tr 



{Pk'ri'^) 



(32) 



where "'^ — ^ q ^ij |0) pl^^ (minus) 

one eigenvector of cFz. 



Figures and |7| plot the bipartite components of 
Eq. (|29|l for the various two-qubit marginals obtained 
after tracing over any one of the three subsystems {A, 
I, or II) as functions of the parameter r, while Fig. |S| 
does the same for the relevant single particle information 
measures. The complementary nature of these quanti- 
ties is illustrated by the pattern/color scheme chosen for 
the curves, where each pattern (color) corresponds to a 
unique set of properties satisfying Eq. . Adding to- 
gether all of the curves that contain a common attribute 
in Figs. El -|H1 yields the constant value one, independent 
of r. For example, taking the sum of the two dashed 




FIG. 6: (Color online) Bipartite tangles as a function of r. 
Thick solid curve (Blue): between Alice and Rob in region 
I. Thin dashed curve (Red): between the modes in regions I 
and II. Dotted curve (Green): between Alice and anti-Rob in 
region II. 




FIG. 7: (Color online) Separable uncertainties as a function of 
r. Thick solid curve (Blue): between Alice and Rob in region 
I. Thin dashed curve (Red): between the modes in regions I 
and II. Dotted curve (Green): between Alice and anti-Rob in 
region II. 

(red) curves in Figs. El and [3 and the two curves contain- 
ing dashes (red) in Fig. |S| thick-thin dashed (blue-red) 
and dot-dashed (green-red), corresponds to the equality 

Tiji + VI Ji + ^{Pi) + ^{pii) = 1- Note that the 
curves in Fig. [HI representing information encoded in the 
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FIG. 8: (Color online) Single qubit properties as a function of 
r. Dotted-solid line (Blue-Green): Alice (always zero) . Thick- 
thin dashed curve (Blue- Red): Rob in region I. Dot-dashed 
curve (Green- Red): anti-Rob in region II. 



individual qubits, each contribute to exactly two distinct 
complementarity relations. 

The tradeoffs expressed by Eq. (|29|l . and illustrated 
in Figs. El - El show how Rob's acceleration determines 
the distribution of individual and bipartite properties 
of the system. For example, adding together the thick 
solid (blue) and dotted (green) curves in Fig. shows 
that, although the entanglement in any bipartite subsys- 
tem changes with acceleration, the total entanglement 
between Alice and (the modes in regions I and II) is al- 
ways maximal and constant. Further, since Rob is forced 
to trace over the causally disconnected modes in region 
II, the entanglement arising between the modes in region 
I and region II is ultimately responsible for the Unruh 
radiation that he sees. Indeed Alice, who has access to 
the modes in both regions, always sees an undisturbed 
vacuum. However, our analysis shows that the maxi- 
mum amount of thermalization that Rob can experience 
(even at infinite acceleration) is limited by the amount 
of entanglement that exists between Alice and Rob in an 
inertial frame. 

One straightforward way to see this is to consider the 
case in which Alice and Rob initially share no entan- 
glement (e.g. an arbitrary product state between Alice 
and Rob). Then, as Rob's acceleration approaches infin- 
ity, we find that tiji approaches its maximum possible 
value of one (see Eq. JTHl) with r = 7r/4 |23|). This is 
in contrast to the situation considered here in which the 
initial entanglement between Alice and Rob is maximal, 
constraining r/j/ to a maximum value of 1/4. 

Since the overall state of our system is always pure, we 
may also apply the result 



-frW-f 52 (pfc) = l. 



(33) 



which holds for an arbitrary pure state of three qubits (a, 
6, and c). The quantity r'*'^ = X^j^^fe T?',^ ^ measure of 
the total pairwise entanglement shared between qubit k 
and all other qubits. Equation (|33|) quantifies an explicit 



tradeoff between the inherent three-body correlations in 
a tripartite system, the total pairwise entanglement of a 
selected qubit, and the single particle properties of the 
qubit. 

Because the three tangle Eq. (|28|l is always zero for 
our system, this pure state complementarity relationship 
simplifies to 



r(fe) 



52 {pk) = 1. 



(34) 



It is then straightforward to see that this tradeoff is also 
captured by the pattern/color scheme used in the above 
figures. Simply adding the sum of any two curves in 
Fig. El to the curve in Fig. |S1 that is composed of the same 
two attributes (or colors) again yields the constant value 
one, regardless of Rob's acceleration. In particular, this 
shows that the entanglement between Alice and Rob in 
region I cannot vanish at infinite acceleration as it did in 
the bosonic case. This is because the entanglement tiji 
generated between regions I and II (dashed/red curve 
in Fig. EJ and the single particle properties S"^ (pi) that 
manifest for Rob (dashed/red-blue curve in Fig. jSJ are 
insufficient to satisfy Eq. H34|l . Instead, the contribu- 
tion from the entanglement / between Alice and Rob 
(solid/blue curve in Fig. EJ must also be taken into ac- 
count to ensure that ta,i + tjjj + S"^ (pi) = 1. 

Finally, we note that the various two-qubit marginal 
density operators for our system have a very specific 
form; they are all examples of the maximally entangled 
states with fixed marginal mixednesses (MEMMS) iden- 
tified in [23. Indeed, this fact is closely related to the 
absence of tripartite correlations in our system, which 
implies that any entanglement that exists in the sys- 
tem must necessarily be bipartite. As shown in 0, the 
MEMMS are characterized by the relationship M (p) = 
rj (p), i.e., the marginal mixednesses are equal to the sep- 
arable uncertainties. Thus, our fermionic system encodes 
the maximum amount of bipartite entanglement consis- 
tent with the separable uncertainty in each two-qubit 
marginal for every possible value of the acceleration rate. 



VII. SUMMARY AND FUTURE DIRECTIONS 

We have studied the behavior of the entanglement be- 
tween two modes of a free Dirac field in a non-inertial 
frame in flat spacetime from the point of view of two 
observers, Alice and Rob, in relative uniform accelera- 
tion. Our results show that entanglement existing be- 
tween Alice and Rob in an inertial frame is progressively 
degraded by the Unruh effect as Rob's rate of acceleration 
increases. However, unlike the bosonic case in which the 
inertial entanglement vanishes in the limit of infinite ac- 
celeration, in this case the entanglement achieves a min- 
imum value of I/V2 for the entanglement of formation 
(and of 1/2 for the tangle). This fundamental difference, 
a consequence of the fact that fermions have access to 
only two quantum levels vs. the infinite ladder of exci- 
tations available to bosons, means that in this case Alice 
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and Rob always share some entanglement which can in 
principle be used as a resource for performing certain 
quantum information processing tasks. Further analysis 
shows that the total (quantum plus classical) correlations 
in the system, as quantified by the mutual information, 
behaves in a manner reminiscent of the bosonic case, de- 
creasing from 2 for inertial observers to unity in the case 
of infinite acceleration. 

Considering the (causally inaccessible to Rob) modes 
in region II to be a third subsystem allows us to analyze 
this system in terms of entanglement sharing. In doing 
so, we find that the overall tripartite pure state never en- 
codes any inherently three-body correlations, regardless 
of the rate of acceleration. Any entanglement existing 
in the fermionic system is therefore necessarily bipartite. 
Such entanglement is known to be an invariant quantity 
for inertial observers, although different inertial observers 
may see these quantum correlations distributed among 
multiple degrees of freedom in different ways. However, 
in this analysis we find no indication that the entangle- 
ment in any bipartite subsystem is conserved when Rob 
is allowed to accelerate. In fact, our results show that 
the presence of the communication horizon, which iso- 



lates the modes in region II from an accelerated observer 
in region I, plays a key role in degrading the inertial en- 
tanglement between Alice and Rob. 

Further insight into this behavior is gained by ap- 
plying both pure and mixed state complementarity re- 
lations to different divisions of the system into subsys- 
tems. The ability of this formalism to identify the dif- 
ferent types of information encoded in a quantum state 
facilitates the study of how various subsystem proper- 
ties depend on Rob's rate of acceleration. For instance, 
the constraints imposed by these relations illustrate how 
multi-qubit complementarity prevents the entanglement 
between Alice and Rob from vanishing at infinite accel- 
eration as it does in the bosonic case. Additionally, we 
find that the amount of vacuum thermalization (due to 
Unruh radiation) that Rob experiences at infinite acceler- 
ation is constrained by the amount of entanglement that 
he shares with Alice in an inertial frame. 

One possible avenue for further research along these 
lines is to study the entanglement between Alice and Rob 
in the case that Alice has acceleration oi and Rob has 
acceleration a2. In this case the density matrix, after 
tracing over region II is 



/ cos^ ri cos^ r2 



PA,I 



\ cos ri cos r2 



9 . 2 

cos ri sm r2 








sin^ ri cos^ r2 



cos ri cos r2 



1 + sin^ ri sin^ r2 / 



where tanri — exp(— wc/oi) and tanr2 — exp (— 010/02). 
In this case the logarithmic negativity is = log2(l + 
cos^ ri cos^ r2). Thus, we see that the entanglement is 
further degraded by having two accelerated observers, 
but again remains finite, in this case taking on the value 
N = log2(5/4) in the infinite acceleration limit. 

Another aspect of this problem that deserves further 
consideration is the nature of the communication that is 
and is not allowed between observers in different regions. 
Of course, no communication is possible between the left 
and right Rindler wedges (regions I and II). However, 
even when two parties are not causally disconnected from 
one another (e.g. Alice and Rob in Fig.QJ, the presence 



of a horizon due to the acceleration of an observer im- 
poses a unidirectionality on any classical communication 
occurring after the inertial observer crosses the horizon. 
Given the connection between the horizons seen by accel- 
erated observers, and the event horizon of a black hole, 
this in turn suggests the potential for gaining further in- 
sight into questions related to the black hole information 
paradox. 
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